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A particle in a 3-D box 



In an ideal gas, we assume that there is little interaction 
between atoms. We shall learn how to calculate the total 
energy and other properties, and apply this to electrons, 
phonons, photons, liquid helium, etc. 

From the kinetic theory of gases, we know that the average 
energy of each atom in an ideal gas is 

1 2 3, 

-mv z = —kuT. 

2 2 B 

This means that the total energy is 

U = -Nk B T 
2 " 

and the heat capacity is 

C = ™ = —Nkfi. 
dT 2 
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Density of States 



A particle in a 3-D box 



We first determine the energy states for a particle in a 3-D box. 

yi 





We solve Schrodinger's equation for the particle: 

h 2 (d 2 ip + d 2 ip + d 2 ip\ +y ^_ E ^ 

2m \dx 2 dy 2 dz 2 J 
where V is zero in the box, and infinite at the walls. So the 
wave function must have zero amplitude there. 
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Energy states of the particle 



This is done in the standard way by trying the solution 

ip = sin k x x sin k y y sin k z z 

where (k x ,k y ,k z ) are unknown constants. They can be 
determined from the boundary condition that the wave function 
must have zero amplitude at the wall. 

In the x direction, suppose the walls are at x = 0 and x = a. 
The condition can be written as: 

sin k x a = 0. 

This gives 

k x a = %7r or k x = , 

a 

where n x is a positive integer. This means that k x is quantised. 
Similarly, for the y and the z directions, we find 

k y = and k z = . 

a a 
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Density of states. 



To find the total energy of N particles, we need to add up all 
their energies: 

U = Y, n i e i- 

where e$ is the energy of level i and is the number of 
particles at that level. 

At room temperature, many of these levels could be occupied. 
This is because for an ideal gas, spacing between energy levels 
is much smaller than the average energy of a particle. 

As a result, a plot of n.; against may look like this. 



Energy states of the particle 



We can then substitute the solution 

ijj = sin k x xs\r\ k y y sin k z z 

into the Schrodinger's equation 

h 2 fd 2 ip + d 2 ip + d 2 ip\ +y _ E ^ 
2m \dx 2 dy 2 dz 2 ) 
to find the energy. This gives 

h 2 
2m ' 

It can also be written in the form 

h 2 k 2 



-{k 2 x + k 2 + k 2 ) = E. 



where 



E 
2m 



k — m 



k is called the wavevector. 



Density of states. 



This suggests that we may approximate the graph to a curve, 
and the sum to an integral. 




Within a small interval de, the energy e is nearly constant. If we 
know the number of particles in this interval, we can multiply 
by e to find the total. 

Because of the quantisation condition, there can only be a 
certain number of states in de. For each state, there is a 
certain probability that it would be occupied. Lets start by find 
the number of states. 
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Counting states 



The wavevectors 

_n x n 
a 

are quantised at uniform intervals in all three directions. 

we can imagine a k space in which the the x coordinate is k x , 
and so on. If we use a point to represent each state, we would 
get a lattice like this. 




Because k x , k y and k z are all positive, only the points within 
1/8 of the sphere of radius k, where k is the wavevector for e. 

Since the spacing between points is n/a, the volume associated 
with each point is (n/a) 3 . Therefore, the total number of 
states is one eighth of the sphere volume divided by (ir/a) 3 : 
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Counting states 



First, we find the number of all states below energy e. This 
energy is related to the wavevector k by 

h 2 k 2 



2m 



k is in turn related to the components by 



k — m 

This describes a sphere in k space. 
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DOS 



We find that the total number of states below k is 

G(fc) = 

where V is the volume a 3 of the box. 

Suppose k is increased by dk. Then the number of states 
increases by dG. So the density of states is 

dG(fc) Vk 2 

9(k) = ^r=2^- 
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DOS 



We have derived g{k) by first deriving G(fe). We can also derive 
gr(fc) directly. 




Consider a spherical shell in k space with a thickness dk. The 
states in the shell lie between wavevectors k and k + dk. 

Each state is associated with volume (n/a) 3 . So the number of 
states in the shell is volume of the shell divided by (n/a) 3 . 
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DOS 



We have earlier defined the DOS in terms of energy, e. 

To convert the variable from k to e, we must use the method 
for probability density function: 

g 6 (e)d£ = g k (k)dk. 

The subscripts are added here to emphasise that g{e) and g{k) 
are different functions. 
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DOS 



Volume of the spherical shell is 4nk 2 dk. We only take 1/8 of 
this because only positive components of k x , k y , k z are needed. 

1 9 / 7T \ ~ 

number of states in dk = — x Auk dk -f- 



a I 



, . number of states in dk 
9{k) = dk 

1 9 /7T\ 3 

q(k) = - X 4nk - 
yW 8 \a 



, H Vk 2 

27T Z 



where volume V = a 3 . 
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Changing variables 

Rearranging, 

de 

We can now substitute the relation 

h 2 k 2 



2m 

to find g £ (e). The answer is 

5 (e) = - ft ^(2me) 1 / 2 

where we have dropped the subscript again, as is the normal 
practice. For comparison: 

, H Vk 2 

27T Z 

If you substitute e for k in g(k), you DO NOT get g(e). g(k) 
and g(e) are DIFFERENT functions! 
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Changing variables 



It is also useful to have the formula for total number states 



in terms of energy. 



Since no e or k intervals are involved, we can simply substitute 

h 2 k 2 

£ 



2m 

"he answer is 

AttV q /i 

G(s) = 3 ^(2me) 3 / 2 
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Ideal gas energy states 



We first show that it is very unlikely for two atoms to occupy 
the same energy state. 

At a temperature T, we know from simple kinetic theory of 
ideal gas that the average energy of a gas atom is about 
3k B T/2. 

First, we need to know how many energy states there are 
below this energy at room temperature. We can use the 
formula we have derived: 

by setting e = 3k B T/2. 
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Most Likely Macrostate 



17 

Ideal gas energy states 



4nV Q /r-Y 

G(3k B T/2) = ^ (2me) 3 / 2 

Putting in the numbers for one mole of ideal gas at room 
temperature, we find that the number of states below 3k B T/2 
is about 10 30 . 

1 mole of gas contains about 10 24 atoms. 

This means we have about 10 30 10 24 = 10 6 states per atom. 
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Ideal gas energy states 



This means we have about 1 million energy states for every 
atom. 



So it is extremely unlikely that two atoms would ever occupy 
the same energy state. 

We are now going to derive the energy distribution of these 
atoms. 
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Energy levels as partitions 



Consider the i th bundle. There are energy states and n, 
atoms. 



Recall that it is extremely unlikely for two atoms to be in the 
same state. 



Think of each state as a partition. The partition and the atom 
are like two different objects. We need to find the number of 
ways to arrange atoms and partitions in a row. 
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Grouping energy levels into bundles 



We shall borrow the idea from histograms - divide the energy 
levels into bins, or bundles. 



g M levels in the bin 
g j levels in the I th bin 
gr,., levels in the (/-1) ,h bin 



Suppose each bundle covers an energy interval de. 



We want to find out how many ways we can arrange a number 
of atoms in a particular bundle. 
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Number of arrangements 



There are m atoms and partitions. 

We can use the combination formula n C r = ,, w! .. . 

rl{n—r)\ 

This gives, for one bundle, ^"'1"„ 9 P ! ways. 

To find the total number of arrangements for all bundles, we 
must multiply together the answer for every bundle: 

i-r ( n i + Qi) ! , - 



Since the number of energy states in each bundle is much 
larger than the number of atoms, >> n it we can simplify this. 
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Expand 



n _ tt (JH + nj)...(9i + l)(gi)-(l) = tt (ffi + n i)---(di + !) 



Since gr; >> rij, (si + n ; ) « 5i and + 1) « Therefore 



This gives the number of microstates in the macrostate 
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Most probable macrostate for the ideal gas 



We apply the Lagrange multiplier method as follows: 

1. Instead of Q, maximise InQ. 

2. Apply Stirling's approximation: 

Inn = Y\n i \r\g i -n i \nn i + n i ] = £ In^/rij)] + AT 

3. The number of particles is fixed: 

iV = J>> 

i 

4. The total energy is fixed: 

U = Yl n i e i- 

i 

5. Define the Lagrange function as In Q + AiAT + X2U. 

6. To find the maximum, differentiate with respect to nf 

In ft - Inn; + Ai + A 2 £i = 0 
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If a macrostate is specified by (n±, 712, n-3, ...), or in short, 
then the above expression gives the number of microstates in 
the macrostate. 

We are interested in the most likely distribution of 
(ni, n 2 , n 3 , ...), or the most probable macrostate. 

We can once again apply the Lagrange multiplier to find this. 
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Most probable macrostate for the ideal gas 
We have found 

Ing; - Inn; + A x + A 2 £i = 0 
Rearranging, we get 

ni = Agiexp(\ 2 £i), 

where 

A = exp(A;L). 

In the same way as before, by using our knowledge of entropy, 
we can show that 

A 2 = - 

k B T 

This gives the Boltzmann distribution again: 

n, = Ag t exp(-ej/fc B T), 
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Approximate to continuous energy 



We are mainly interested in the energy range of the order of 
k B T (which means up to a few times of that). Most of the 
atoms would fall in this range. 

We have seen previously that there is a very large number of 
energy states in this range. 

Since the energy levels are so closely spaced, we can 
approximate them to continuous energy. 

We have previously obtained the density of states, which is the 
number of states per unit energy, g{e). 

Consider the i th bundle. Assume that this occupies an energy 
interval of de. 

The number of energy states is then given by eg = g(e)de. 
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Approximate to continuous energy 



Substituting the formula for g(e), we get 

n(e) = A Am ^ 3 V (2me) 1 / 2 exp(- Si /k B T) 

Note that the total number of atoms is then given by 
integrating over all energy: 

roo 

N = / n(e)de 
Jo 

This relation is known as the normalisation condition. 
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Approximate to continuous energy 



In a similar way, we can define a number density n(e) as the 
number of particles per unit energy. 

So the number of atoms in the bundle would be 

nj = n(e)de. 

Therefore, the distribution that we have obtained, 

rii = A gi exp(~/3€i), 

can be rewritten as 

n{e)de = Ag(e) exp(-£i/k B T)de 
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Maxwell-Boltzmann distribution 



The unknown constant A can be determined by integrating the 
normalisation condition (with the help of the table of integrals): 

3/2 



V \2nmk B TJ 



This gives 



n ( £ ) = 7 , 27 ^3/ 2 ( £ ) 1/2 ^xP(-e»/fcffT) 



The energy is just the kinetic energy, e = mv 2 /2. So 

/2\ 1/2 / m \ 3 / 2 
n (v) = N [iT^p) exp(-mw 2 /2k B T)v' 

This is the Maxwell-Boltzmann distribution of speeds. 
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When would it break down? 



We have assumed that: "It is extremely unlikely that two 
atoms would ever occupy the same energy state." 



Recall that we have arrived at this by looking at the number of 
states below the mean energy 3fc^T/2, when T is at room 
temperature. 

If T is very small, there may be far fewer states below the mean 
energy. Then it would be very likely for two atoms to occupy 
the same state, and the Maxwell-Boltzmann distribution would 
not be valid. 

Instead, quantum statistics - such as Fermi-Dirac or 
Bose-Einstein - have to be used. We shall learn about these 

later. 
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Boltzmann constant 



We can now return to the earlier lectures on derivation of the 
Boltzmann distribution 

rii = Aexp(-ei/k B T) 
and tie up a loose end. 

We can now find the total energy in an ideal gas. Replace the 
sum 

U = Yl n i £ i 

by the integral 

roo 

U = n(e)ede. 
JO 

Substitute the formula for number density: 

roo 2irN 1 /o 

U= ( £ ) 1 /2 e xp(^ £ /fc J? r) £ d £ . 

JO (■KkaT) S ' 2 
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Boltzmann constant 
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Boltzmann constant 



Integrating this gives the formula for the total energy that we 
know from kinetic theory of gases: 

3 

U = -Nk B T 
2 u 

This is evidence that the Boltzmann distribution for the ideal 
gas is correct. 



Differentiating with respect to T gives the heat capacity: 



The measured value is: 



C = -Nk B . 
2 a 



3 

C = -nR. 
2 



R is the gas conant and has the value of 8.314 J mol 1 K 



i k—i 
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Recall that k B is related to one of the Lagrange multipliers. 
We can find it by comparing the measured C above, with the 
formula derived using the Boltzmann distribution. 

Comparing 

C = -Nk B 

and 




we find 

Nk B = nR. 

Since the number of moles n is related to the number of 
particles N by 

N 



Worked Examples 
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we have 

N 

Nk B = R 

N A 

and finally, we obtain the value for k B : 

t R 

N A 

which is indeed the Boltzmann constant as we have assumed. 



Example 1 



Show that the partition function of an ideal gas: 

Z= ^—(2m£) 1 / z exp(-£/k B T)dE 

JO 

is given by: 

'2-Kmk B T\ 3 / 2 



Z = V 



h 2 



[You are given that 

° o 2 2 

n\(2a) 2n a 



r ^e-Vd* - (2n)! " 1/2 

J—OO 
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First, move the constant factors outside the integral: 

Z = ^— (2m) 1 / 2 / £ 1 ' 2 exp(-e/k B T)de 



h 3 

There is no power of 1/2 in the integration formula that is 
given. So we start with the substitution 

£ = X 2 

Then 

„ 4mirV 



(2m) 1 / 2 f°° xexp(-x 2 /k B T).2xdx 
JO 

V (2m) 1 / 2 2 [°° x 2 exp(-x 2 /k B T)dx 



h 3 



h 3 Jo 
Next, set n to 1 in the given formula 



oo _ 2 2 27T 1 / 2 

x 2 e~ a * dx = ^Lg- 



The lower of limit of this integral is — oo. To make this agree 

with the integral for the partition function Z ... 
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We get 

Z = -^(2m) 1 / 2 2— (k B T) 3 ' 2 

Simplifying, 

z = v( 2m l^ T f 2 

we get the formula for the partition function. 
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We change the lower limit to 0 and halve the right hand side, 
and the formula becomes: 



oo 2 2 

x 2 e~ a x dx 



7T 1 / 2 



(2a) 2 a 

Now comparing with the partition function again: 

Z = 4m !" y (2m) 1 / 2 2 [°° x 2 exp(-x 2 /k B T)dx 
h 6 JO 

the integrals become the same if we set 

2 1 

a = . 

k B T 

Substituting this into the formula gives 

rV2 



J"x 2 e- 2 ^ T dx = ^—{k B T) 3 l 2 
Substituting this into partition function intergral, 
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Example 2 

What is the root-mean-square speed of a helium atom at room 
temperature? 

[Given that the relative atomic mass of helium is 4. 
Atomic mass unit u is 1.6605 x 10~ 27 kg. 
Boltzmann's constant k B is 1.3807 x 10" 23 J K -1 . ] 
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First, we find the mean of v 2 . This is can be obtained by 

v n\e)de. 



i r°° 2 



e = —mv 2 . 



N Jo 

We have multiplied v 2 by the number particles at each energy 
interval de, and divided by the total number N. 

The energy of the atom in an ideal gas is just the kinetic 
energy, 

1 

—j 
2 

If we rewrite the above integral for v 2 in terms of the energy, 
we get 

- 2 f°° 

v 2 = / en(e)d£. 

NmJo 

This integral is just the total energy U, i.e. 

v 2 = —U. 
Nm 
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We can now calculate the root-mean-square speed. The mass 
of the helium atom is 

Au = 4 x 1.6605 x 10" 27 = 6.642 x 10" 27 kg. 
Taking the room temperature as 298 K, the answer is 

V m 

3 x 1.3807 x 10" 23 x 298 
~ \ 6.642 x 10" 27 

= 1363 m/s 
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We have previously obtained the energy U from the partition 
function: 

3 

U = -Nk B T. 
2 a 

So the mean square speed is 

o 2 3 AT , m 3k B T 

v 2 = Nk B T = — — 

Nml m 

The root-mean-square speed is then 



3k B T 
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